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Abstract. A vanishing viscosity method is formulated for two-dimensional transonic 
steady irrotational compressible fluid flows with adiabatic constant 7 £ [1,3). This 
formulation allows a family of invariant regions in the phase plane for the corresponding 
viscous problem, which implies an upper bound uniformly away from cavitation for the 
viscous approximate velocity fields. Mathematical entropy pairs are constructed through 
the Loewner-Morawetz relation by entropy generators governed by a generalized Tricomi 
equation of mixed elliptic-hyperbolic type, and the corresponding entropy dissipation 
measures are analyzed so that the viscous approximate solutions satisfy the compensated 
compactness framework. Then the method of compensated compactness is applied to 
show that a sequence of solutions to the artificial viscous problem, staying uniformly 
away from stagnation, converges to an entropy solution of the inviscid transonic flow 
problem. 



1. Introduction 

In two significant papers written a decade apart, Morawetz [271 I28| presented a program 
for proving the existence of weak solutions to the equations governing two-dimensional 
steady irrotational inviscid compressible flow in a channel or exterior to an airfoil. As is 
well known, the classical results of Shiffman and Bers [I] apply when the upstream 
speed is sufficiently small, for which the flow remains subsonic and the governing equations 
are elliptic (also see [T^l EH ESI IIHl EZl HH1 IIHl HOI)- However, beyond a certain speed at 
infinity (determined by the flow geometry), the flow becomes transonic which, coupled 
with nonlinearity, yields shock formation (cf. |26|). Morawetz's program in |27l I28| was 
to imbed the problem within an assumed viscous framework for which the compensated 
compactness framework (see Section 6 of this paper) would be satisfied. Under this as- 
sumption, Morawetz proved that solutions of the as yet unidentified viscous problem have 
a convergent subsequence whose limit is a solution of the transonic flow problem. 

The purpose of this paper is to present such a viscous formulation, hence completing 
part but not all of Morawetz's program. Specifically, a vanishing viscosity method is for- 
mulated for two-dimensional transonic steady irrotational compressible fluid flows with 
adiabatic constant 7 6 [1,3) to ensure a family of invariant regions for the corresponding 
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viscous problem, which implies an upper bound uniformly away from cavitation for the vis- 
cous approximate velocity fields. Mathematical entropy pairs are constructed through the 
Loewner-Morawetz relation by entropy generators governed by a generalized Tricomi equa- 
tion of mixed elliptic-hyperbolic type, and the corresponding entropy dissipation measures 
are analyzed so that the viscous approximate solutions satisfy the compensated compact- 
ness framework. Then the method of compensated compactness is applied to show that 
a sequence of solutions to the viscous problem, staying uniformly away from stagnation, 
converges to an entropy solution of the inviscid transonic flow problem. 

On the other hand, Morawetz's assumption of no stagnation points for the flow has not 
been removed; and the slip boundary condition (u, v) n = on the obstacle is only satisfied 
as the inequality (u, v) ■ n > for the general case (the desired condition (u, v) ■ n = 
may be achieved when a solution has no jump with certain regularity along the obstacle), 
where (u, v) is the fluid velocity field and n is the unit normal on the obstacle pointing 
into the flow. Both issues may reflect rather complicated boundary layer behavior for 
transonic flow [311] and require further investigation, which are the topics for subsequent 
research. In addition, we note that, when 7 > 3, cavitation is indeed possible and will be 
the topic of a sequel to this paper. 

This paper is divided into nine sections after the introduction. In Section 2, the classical 
fluid equations are presented for isentropic and isothermal irrotational planar flow. In 
Section 3, the fluid equations are rewritten in the polar variables (u,v) = (q cos 6, qsin9). 
In Section 4, we analyzes the behavior of the Riemann invariants in the supersonic region, 
which guides the design of artificial viscous terms in our vanishing viscosity method. In 
Section 5, we continue our analysis and set boundary conditions for the viscous system. 
This yields a family of invariant regions in the (u, v) fluid phase plane for 7 £ [1, 3), where 
7 is the ratio of specific heats if 7 > 1 and 7 = 1 denotes the case of constant temperature. 
In particular, any invariant region is uniformly bounded away from the cavitation circle 
in the phase plane, which yields a upper bound uniformly away from cavitation in the 
viscous approximate solutions. In Section 6, we formulate a compensated compactness 
framework for steady flow. In Section 7, we first construct all mathematical entropy 
pairs for the potential flow system through the Loewner-Morawetz relation by entropy 
generators governed by the generalized Tricomi equation of mixed elliptic-hyperbolic type. 
Then we introduce the notion of entropy solutions through an entropy pair by a convex 
entropy generator suggested by the work of Osher-Hafez- Whitlow |32) , In Section 8, we use 
the previously derived uniform L°° bounds on (u e ,v e ), plus the convex entropy generator 
introduced in Section 7 to guarantee our problem is within the compensated compactness 
framework. In Section 9, we develop Morawetz's argument in [2S] for proving convergence 
of a subsequence of solutions of our viscous problem to the irrotational fluid problem. 
Finally, in Section 10, we prove the existence of smooth solutions to our viscous problem. 



2. Mathematical Equations 
Consider the artificial viscous system in a domain del 2 : 
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where R\ and R 2 are the artificial viscosity terms to be determined, and (u, v) is the flow 
velocity field. For a polytropic gas with the adiabatic exponent 7 > 1, p = p(p) = p 7 /7 is 
the normalized pressure, the renormalized density p is given by Bernoulli's law: 

p = p(q) = {l-^-±q 2 )^, (2.2) 

where q is the flow speed defined by q 2 = u 2 + v 2 . The sound speed c is defined as 

c 2 =p'(p) = l-^q 2 . (2.3) 

At the cavitation point p = 0, 

Q Ql,ll l \ ! > ■ 

7 - 1 

At the stagnation point q = 0, the density reaches its maximum p = 1. Bernoulli's law 
(t2~2l is valid for < q < 

q C av At the sonic point q — c, (|2,3j) implies q 2 — ?^ . Define 

the critical speed q cr as 



7 + 1 

We rewrite Bernoulli's law (12. 2|) in the form 

Thus the flow is subsonic when q < q cr , sonic when q = q cr , and supersonic when q > q cr . 

For the isothermal flow (7 = 1), p = c 2 p where c > is the constant sound speed, and 
the density p is given by Bernoulli's law: 

, u 2 + V 2 s 
p = p exp( 2&~' 

for some constant po > 0. In this case, q cr = c. After scaling, we can take c = po = 1. 

3. Formulation in Polar Coordinate Phase Plane 

We now use the polar coordinates in the phase plane: 

u = qcos9, v = qsin9, 

and rewrite the viscous conservation laws (|2,lj) in terms of (q, 9). Write the second equation 
of JZB as 

p x u + pu x + p y v + pv y = R 2 

or 

p'{q)q x u + pu x + p'{q)q y v + pv y = R 2 , 

and use 

'/ V ''" • ' </.r ~ ' 

to find 



qr = a/u 2 + v 2 , q x = -(uu x + vv x ), q y = -{uu y + vv y ), p'(q) = — - 



c 2 



(c - u )u x - uv(y x + u y ) + (c — v = — ii 2 - 
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Then 1(2.1(1 becomes 

r -Ri 



A 



+ B 



-R 2 



where 



A 



c — u —uv 



D 



-uv c — V 



Thus, in terms of (q,9), we obtain 



where 



sin# 



^!cos( 

c^q 



+ B 1 



-Ri 
m R2 



J cos I 
sin 8 



cos 8 
— j-^- sm( 

c A q 



-qsml 
cos 8 



In terms of (p, 9), using 



Px = PqQx 



PQ 

- 2 Py 



pq 

'2 q y> 



and thus 



we find 



q 




c 2 

pqP x 


5 


q 




c 2 

pqPy 


e 


X 


Ox 




e 


y 





An 



+ B 2 



-Ri 
-R 2 



where 



sin ft — qcost 

1 (c 2 — g 2 ) cos ft pq sin£ 



Bo 



±(c 2 -g 2 )sini 



We symmetrize 1)3.2(1 to obtain 



-q sin ft 
-pq cos ( 



A, 



+ #3 



-Ri 

q 2 

^iR 2 



where 
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4. Choice of Artificial Viscosity from an Analysis of Riemann Invariants 

The two matrices A\,B\ in H3.1|) commute, thus their transposes commute and they 
have common eigenvectors. The eigenvalues of A\ and B\ are 



X± = — sin 6 ± — cos 9, fi± = cos ± sin ( 

c c 



The left eigenvectors of A\ are 



and thus the Riemann invariants W± satisfy 



8w± <w ± y^^2 

-55- = 1, = =F for q > c. (4.1 

OP ag gc 

Multiply (EH) by (^, ^) to obtain 

dw± dw± 0W± , 1 aW± 

A±^ h p±^ — = tt~ R\ H (4.2) 

ox ay oq pq ot> 

We now consider the viscosity terms of the form: 

Rx = eV ■ (ai(p, 0)V0) , i? 2 = eV • (a 2 (p, 0)Vp) . (4.3) 

In particular, for a special choice of o\ and 02, we have the desired relations for the 
Riemann invariants. 



Proposition 4.1. // 



c 2 

o-l = l, o 2 = 1 k for q> c, (4.4) 



i/ien the Riemann invariants W± satisfy the following equations: 

qc f. dW± dW±\ . t t 7 (7-3)g 2 + 4c 2 ,^ l2 fA „. 

for 7 > 1. 

Proof. We first focus on the "+" part, since the "-" part can be done analogously. The 
proof is divided into three steps. 

Step 1: Substitute (pL"T|) into to obtain 

A+ ^± + , + ^ = VZZZ fil+ i fi , (4 . 6) 

ax ay qc pg 

Multiplication of (|4.6j) by gc gives 

\/q 2 -c 2 

qc ( dw+ dw + \ dw +v c 



v /^3^2 v T ax ^ d y J de • py /, 

^ = —pq/c 2 , we have 

(9W+ dW + dq c^J q 2 — c 2 

op og dp pg 2 
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and then 

c q 2 dW + 



Therefore, 



pyV - c 2 q 2 ~c 2 dp 



- c 2 V dx ^ dy ) dO q 2 c 2 dp 

For the choice of viscosity terms in (|4.Hj) , we obtain 
qc / dW+ dW^ 



V + dx ^ dy 

= e^V ■ (a^p, 0)V0) + e -^^ V • (a 2 (p, 0)Vp) 
dU+ ( ' do\„ _„ 9o"i,„„,o\ , „ N cW+ .„ 



(4.8) 



g 2 0W+ /da 2| ^ l2 3ct 2 ^ „„\ , „, q 2 dW+ 



.Step <?: From a direct calculation, we know 

A^ = ^|Vp|* + ^Ap + A0, 

then substitution for A# in ()4.8j) gives 
qc / 

A+— h/U H 



' — iVp • V0 + -^|V0| 2 + ecx^p, 0)— t - -^Vp) 2 - -^Ap 



dO \dp ' dO 1 1 / lvr ' 7 90 V ^p 2 1 ri dp 

g 2 dW+ (da 2 ,„ l2 dcio —A , „, g 2 dW+ 

To eliminate Ap, we choose 

„2_ c 2 



0"2 (p) = 0"! (p) 2 

where o\ and cr 2 are independent of 0. Then we have 
qc ( . dW+ dW + 



7^23^2 V ^ dx ^ dy 
dai 9W+ V X70 + £a 

dp dO ' dO \ dp 2 1 " ri J dp q 2 — c 2 dp 



So-icW. ^ n M 3W + / AII , d 2 W +l ^ l2 \ cV 2 g 2 dTy +l ^ l2 
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Finally, write V0 = VW+ - so that 

qc ( . dW+ dW + \ da 1 dW+^ „„, dW+ Ar „ 



^2 _ c 2 \ ^ dx ^ dy J dp 89 r 50 
A convenient choice of ci and 02 is as in (|4.4|l : 

g 2 - c 2 
g 2 

Thus, using = 1, we obtain 

gc /, dW+ dW + \ . /d 2 PF+ d(7 2 g 2 <9W+ 



yS q 2 - c 2 ^ + & r+ 9i/ ^ + V dp 2 dp q 2 - c 2 dp 

Step 3: We now compute the term 

d 2 W+ da 2 q 2 dW + 



dp 2 dp q 2 — c 2 dp 

in gllUJ). 

First we check the simple isothermal case (7 = 1) for which 

p = e 2 , c = 1; — = -pg. 

ag 



Then 



dw+ _ ^ 

dq 



and 



0"2 



_ cM^dg _ yV - 1 
dp dq dp pq 2 ' 
d 2 VF + _ -g 4 + 2g 2 - 2 
dp 2 ~ p 2 q*^ 2 ~^l ' 

7 2 — 1 1 do"2 do 2 dq 



Thus, 



g 2 ' dp dq dp pq A 

d 2 W+ _ da 2 q 2 dW+ ^ q 2 -2 

~dp^ " ~dp~q 2 -c 2 dp ~ pVTT^T' 
and (|4.1U|) becomes, in the case 7 = 1, c = 1, 

dW + dW + \ q 2 -2 , 



\ T 3i ^ dy J ~ pZq>^p~i 

We now consider the case 7 > 1 for which 

2 _ n 7 - 1 2 _ q 2 -c 2 7 + 1 1 



1 , cr 2 



q 2 
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Then 





8W. 


- 




dq 




8W+ 


dW + dq 


v> 


dp 


dq dp 






do2 dq 


2c 2 


dp 


dq dp 


pq 4 ' 


d 2 W + 


d (dW + \ 


dq 


dp 2 


dq\ dp J 


1 dp 



yV - c 2 dp_ _ _pq 

qc dq c 2 



7 2 — r >2 / ^2 "\ n / q2 — r^2. 



— cr / c \ C\l q" — c 



qc \ pqj pq 2 



„2 



dW, 




Thus. 



pq dq \ dp 

c ((3 - 7)g 4 + 2( 7 - 3)q 2 c 2 + 4c 4 
~ 2p 2 q i ^f^ 

d 2 W+ da 2 q 2 dW+ 
dp 2 dp q 2 — c 2 dp 

c , a . „, „, o o , i , <r <-\/ <r - <- 



((3 - 7 )g 4 + 2( 7 - 3)q 2 c 2 + 4c 4 ) + 



2c 2 q 2 C\J q 2 - " 2 



2p 2 (7 4 v / ^ 2 ^ 2V 7 pqU 2 ~c 2 pq 2 

c ((7 - 3)g 2 + 4c 2 ) 
2p 2 q 2 \J q 2 — c 2 
and (|4.10|) becomes, in the case 7 > 1, 

qc / cW+ cW+\ t t 7" (7 - % 2 + 4c 2 .„ l2 

+ - eA ^+ = - £c » o , %-^ |Vp| 2 - (4.12) 



yftf^c 2 \ dx ^ dy J T 2p 2 q 2 ^q^c 2 

This equality is consistent with the case 7 = 1, and this lemma is proved for W+. The 
computation for W- is done analogously. □ 

5. Boundary Conditions, Invariant Regions, and L°° Bounds 

In this section, we consider the level sets of W± to assign appropriate boundary condi- 
tions and identify a family of invariant regions for the viscous problem Q2.1J1 with (|4.3|) - 

First we discuss the isothermal case: 7 = 1 and c = 1. From Proposition 14.11 we find 
that, when q > y/2, W+(x,y) := W+(q(x,y),9(x,y)) cannot have an interior minimum 
since AW + < 0. Similarly, W-(x,y) := W-(q(x,y),9(x,y)) cannot have an interior 
maximum since AW^_ > 0. If we assume for the moment that W+ cannot have a minimum 
on the boundary of our domain, then a solution of the viscous problem which crosses from 
q < \pl to q > v2 at (xq, yo) must satisfy 



W+(x,y) > ^+(x ,yo)| 



--V2- 



Similarly, if we assume for the moment that W- cannot have a maximum on the boundary 
of our domain, then a solution of the viscous problem which crosses from q < y/2 to q > v2 
at (xo,yo) must satisfy 



W.(x,y) < W-(x ,y )\ 
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From the definition of W±, we have 



dW± y/q 2 - 1 
dq q 

Using the substitution q = sect, we can integrate to find 

W± = 9 T vV ~ 1 ± arccos(g~ 1 ) + C± = T (vV - 1 - l) ± (arccos^ -1 ) - j 



where the constants C± is set to be C± = (±y c/ 2 — 1 =F arccos(g 1 ))| g=v ^ = ±(1 — f )• 
Along the level curves W± = const., we have 

d9 _ ± y / g^ r l 

Thus, on the level set of W+, % > and B is increasing when q is increasing; while, on 
the level set of W-, ^- < and 9 is decreasing when q is increasing. 
If 0(xo, 2/0 ) = 9(^0) 2/0 ) = \/2 and we leave q = y/2, we have 

Ofay) - 9 > (vV - 1 - l) - (arccos(g _1 ) - ^ (stay below W+), 



0(x,y) - Go < — \ \/ q 2 — 1 — lj + ^arccos(g 1 ) - —J (stay above W-), 
that is, 

q 2 — 1 — 1^ — ^arccos(g _1 ) — — J < \0(x,y) — 0q\ inside the "apple" shaped region. 

See Figure n with q cav = oo when 7 = 1. 

The same situation occurs when the level set curves W± = W±(qo, Oo) for 0{xo, yo) = Oo 
and q(xo,yo) = Qo > V%, for which we obtain similar invariant regions of "apple" shape 
past the point (qo,0o) in the (u, t;)-plane. 

Now we return to the issue of boundary conditions. Denote the boundary of the bounded 
domain by <90, the boundary of the obstacle by dQi, and the far field boundary by d&2- 
Thus, d£l = d£l\ U 80,2- Since we do not want W+ to have a minimum on boundary d£l\ 
and W- to have a maximum on boundary dVt\, we require 

dW+ dW^ 

— — < 0, — — > at all boundary points, 
an an 

where n is the unit normal into the flow region on d£l. Recall 

q 2 -l 

VW± = V0t Vg. 



Therefore, if we set 
then 



V6>-n = on dn h 



sign(VW± • n) = =F sign(Vg • n) 

at those boundary points where q > 1, i.e. where the Riemann invariants are defined. 
Hence, one resolution of the boundary condition issue is to set 

eo"2V/9 • n = — \p (u, v) ■ n| on 8Ui, 
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Figure 1. Invariant regions of "apple" shape. 



and 

(u, v) — (uoojVoo) = on <9S7 2 , 

with qoo = \(uoo,Voo)\ S (0, q ca v) = (0, oo). Trivially, W± are not even defined on d&2 and 
hence cannot have a maximum or minimum there. Since 

Vp = Vq^- = -pqVq, 
aq 

we have automatically 

sign(VW / ± • n) = ±sign(V/9 ■ n) = =pl on all boundaries, 

and our minimum principle for W+ and maximum principle for are indeed valid. 
Secondly, they formally yield 

p (u, v) ■ n = on 9f2i, 

and 

(u,v) - (uocVoo) — >■ as x 2 + y 2 ^cx), 

if e — > 0. This is the case when any shock strength is zero at its intersection point with the 
boundary as conjectured for the shock formed in the supersonic bubble near the obstacle 
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(see Morawetz 1291). Finally, we see that the relation between o\ and 0*2 is only necessary 
for q > \/2. Any smooth continuations of o\ and a 2 inside q = \[2 with o\,(J2 > suffices. 

With this, we conclude that (q, 0) stay inside the "apple" region for 7 = 1. More 
generally, for 1 < 7 < 3, we have 

Theorem 5.1. Consider the viscous problem 

\v x -u y = eV ■ (oi(/9)V0) , 



(pu) x + (pv) y = e\7 ■ (cr 2 (p)Vp) , 



(5.1) 



with 

q 2 - c 2 

ai = 1, cr 2 = 5 — for q > V2q cr , 1 < 7 < 3, 

q 2 

and the boundary conditions: 

{V6»-n = on dSlt, 
eo" 2 V/9 ■ n = — \p (u, v) ■ n| on dft\, (5-2) 
(u,v) - (uocVoo) = on dQ 2 with q^ < q cav , 

Then all solutions to (|5.1|) - ()5.2|) are bounded on the domain f2 staying uniformly in e 
away from cavitation, i.e. there exists q* < q cav such that q £ < q* that is equivalent to 
p(q 6 ) > P > for some p = p{q*) > 0. More specifically, when 1 < 7 < 3, the viscous 
approximate solutions stay in a family of "apple" shaped invariant regions as shown in 
Figures 1 and 3-5. 

Proof. The theorem has been proved in the above for the isothermal case 7 = l,c = 1, 
and q cr = 1. 

For the case 1 < 7 < 3, from (|4.5|) . we first require 

(7 - 3)q 2 + 4c 2 < 0, 
which implies q 2 > ^^c 2 . Then, from Bernoulli's law (|2.4ft . we have 



lcr *~ ~C 



3-7 3-7 v 2 

thus q 2 > 2q 2 r , i.e. q > V2q cr . 

We first consider the case that the far field speed qoo < \[2q cr . Then, from the definition 
of the Riemann invariants ()4.1|) . we follow Landau-Lifshitz |22| . page 446, and have 

W± = 9 T (W(q) - W{V2q cr ) 

where 



7 + 1 . /7-l/g 2 



W(q) = \ arcsin * — - — I 1 1 — arcsin 



7-1 V 2 \q, 



P- 

cr 




The two level set curves of W± intersect to form a invariant region of "apple" shape 
past {\p2q cr , 9o) as in Figure Q if 

W{q cav ) - W{V2q cr ) > it, 
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1.5 2 2.5 3 



Figure 2. The graph of 0(7) 

that is, 

a ( 7 ) : = W(q cav ) - W(V2q cr ) 

( h + 1 , \ it 1 h + 1 . /t - 1 . h + l\ 
= [ \ 1 \ / arcsm \ / arcsm \ / > 7r. 

vy 7-1 ' 2 W 7-1 V2 V 4 / 

The graph of the function 0(7) is shown in Figure |21 which shows that a( 7 *) = 7r for 
some 7* ~ 1.224 and 0(7) > n for 1 < 7 < 7*. 

For general 7 S [7*, 3), the level set curves of W± = Oq past (V2q cr ,6o) end at some 
points on the cavitation circle q = q cav and do not intersect each other. We denote 
A(y/2q cr ,9o) the round "apple" shaped region formed by W+ > 9q, W~ < 9q, and the 
cavitation circle q = q cav ', see Figure 3. 

When 0(7) G (f ,ir], we find that, for any #0 G [0, 27r), 

A(v / 2g cr , «o) n A(v / 2g cr , 7T + ff ) 

forms an invariant region for the viscous solutions; this follows from application of our 
minimum and maximum principle for W+ and W— , respectively, at the crossing point on 
the new level set curves. Such invariant regions stay away from the cavitation circle (see 
Figure 4). 

When 0(7) G (f , §], then, for any 6 G [0,2vr), 

A(V2q cr , Q ) n A(V2q cr , | + ) D A(V2q cr , it + Q ) D A(V2q cr , + O ), 

forms an invariant region for the viscous solutions, staying away from cavitation; see Figure 
5. 
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Figure 3. Round "apple" shaped region A(^/2q cr ,9 ) when 0(7) < it 

In general, when 0(7) G (^, 2 n-i )i n = 1)2,... , for each #0 G [0, 27r), it requires an 
intersection of at least In round "apple" shaped regions including A(\^2q cr , 9$) to form an 
invariant region staying away from the cavitation circle. 

When the far field speed (fco G (V2q C r,Qcav), then the level set curves of W± = 
W±(qo,9o), for some qo € {qooilcav) and 9q = arctan^^} , past the point (<?o>#o) either 
intersects each other or end at some points on the cavitation circle q = q ca v As before, we 
denote A(qo, 60) the "apple" shaped region formed by W + > W+(qo, 60), W~ < W~ (qo, 9q), 
and the cavitation circle q = q cav similar to either Figure 1 or Figure 3. Then we can sim- 
ilarly obtain the invariant regions which consist of either a single "apple" shaped region 
or some unions of certain number round "apple" shaped regions including A(q , 6 ), which 
stay away from the cavitation but include the state («ocn^oc) (cf. Figures 5-6). □ 

6. Compensated Compactness Framework for Steady Flow 

Let a sequence of functions w £ (x,y) = (u £ ,v £ )(x,y), defined on open subset C M 2 , 
satisfy the following Set of Conditions (A): 

(A.l) q £ (x,y) = \w £ (x,y)\ < a.e. in Q, for some positive constant q* < q cav < 00; 

(A. 2) d x Qi±(w £ )+d y Q2±(w £ ) are confined in a compact set in H^{Q), for any entropy- 
entropy flux pairs (Qi,Q2) so that (Qi±{w £ ), Q2±{w £ )) are confined in a bounded set 
uniformly in L^ C (Q). 
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Figure 4. Invariant regions A(V2q cr ,6 ) n A(^q cr ,7r + 6 ) when 0(7) G (f ,vr) 

Then, by the Div-Curl Lemma of Tartar |36| and Murat [30| and the Young measure 
representation theorem for a uniformly bounded sequence of functions (cf. Tartar |36| : 
also Ball [!.), we have the following commutation identity: 

< i/(w),Qi + (w)Q 2 -(w) - Qi-(w)Q 2+ (w) > 
=< u(w),Qx + (w) >< v(w),Q 2 -(w) > - < v(w),Qi-(w) >< u(w),Q 2+ (w) >, 

where v = u x>y (w),w = (u,v), is the associated family of Young measures (probability 
measures) for the sequence w £ (x,y) = (u £ , v £ )(x, y). This is equivalent to 

< v(w) <g) !/(«/), I(w, w') >= 0, (6.2) 

where v{w) ® v{w') is a product measure for (w,w') G R 2 x R 2 and 

i(w,w') = (Qi+(«;)-Qi + (^0)(Q2-H-Q2-(^0)-(Q2+H-Q2 + (^0)(Qi-H-Qi-(^))- 

The main point for the compensated compactness framework is to prove that v is 
in fact a Dirac measure by using entropy pairs, which implies the compactness of the 
sequence w 6 (x,y) = (u £ , v £ )(x, y) in L\ oc {VL). Some additional references on compensated 
compactness method include Chen EJ, Dafermos JU], DiPerna jllj . Evans |13j . and 
Serre jHl]. Theorem 15.11 shows that (A.l) is satisfied for our viscous solution sequence 
w £ = (u £ ,v £ ). In the next two sections, we show how (A. 2) can be achieved. 



7. Entropy Pairs via Entropy Generators 
We first construct all mathematical entropy pairs for the potential flow system. 
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Figure 5. Invariant regions A(v / 2<? cr , 0o) n A(v / 2g C r, | + 6 )nA(V2q cr ,ir + 
Q ) n A(V2q CT , f + ) when 0(7) G (f , |] 



For some function V(p, 0) to be determined, multiply (|3.3|) from left by (Vg,V p ) to get 
the entropy equality 

q 2 

Qlx + Q2y = -VgR l + -^-^V p R 2l (7.1) 

c q 

where (Qi,Q2) is defined by 

dQi c 2 aQi pq 3 sing 

— — = — sin Vg - gcos V p , = -g cos 6* 5 5- V p , 

op pq 00 c z — q z 

0Q2 c 2 <9Q 2 q -it . pq 3 cos6 

— — = cos - g sm V p , -—— = -q sin + —5 5- V>. 

op pq 00 & — q z 



(7.2) 



We note that, the appearance of the term c 2 — q 2 in the denominator of l|7.2|) is only 
a consequence of our formalism. It will cancel out as we proceed. Using = 
« = 1, 2, from Q7.2[) . we see that V satisfies the Tricomi type equation of mixed type: 

2 / 3 \ 

— V ee + qV P + ( -P—zV p ) =0. (7.3) 
pq r \c 2 -q 2 F J p 
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Thus, we have defined an entropy pair (Qi,Q2) generated by V. Alternatively, we can 
define an entropy pair (Qi, Q2) generated by H, where H and V are related by 

pH^e -H = -V e , H„ + -H ee = ~^-^V p , (7.4) 

p c z — q A 

with p = p(p) defined by p'(p) = c 2 /q 2 , and H is determined by the generalized Tricomi 
equation: 

H^ + \(l-M 2 )H ee = 0, (7.5) 



and M = q/c is the Mach number. 

Lemma 7.1. The entropy pairs (Qi,Q2) are given by the Loewner-Morawetz relation: 

Qi = pqH^ cos 9 — qHg sin 9, Q2 = pqH^ sin 9 + qHg cos 9, (7-6) 
where the generators H are all solutions of ()7.5|) . 

This can be seen by differentiation of (|7.6|) with respect to (p, 9) and comparison with 

C2J. 

A prototype of the generators H, as suggested in [21], is 

#* = y + y J ^(M 2 -l)dp'dp, (7.7) 

which is a trivial solution of the generalized Tricomi equation (|7.5j) . Notice that H* is 
strictly convex in (p, 9) in the supersonic region. We denote by (Q*, Q%) the corresponding 
entropy pair generated by the convex generator H* . With this, we introduce the notion 
of entropy solutions. 

Definition 7.1 (Notion of Entropy Solutions). A bounded, measurable vector function 
w(x,y) := (u,v)(x,y),q(x,y) < q ca v, is called an entropy solution of the potential flow 
system in a domain £1: 

{V x -Uy=0, 
\(PU) X + (pv)y = 0, 

if w(x, y) satisfies (|7.8|) and the following entropy inequality: 

Q* lx + Q* 2y <0 (7.9) 

in the sense of distributions in $7, in addition to the corresponding boundary conditions 
in the trace or asymptotic sense on d£l. 

The physical correctness of the entropy inequality (|7.9j) is provided by Theorem 2.1 of 
Osher-Hafez- Whitlow [S2J- 

8. H~ l Compactness of Entropy Dissipation Measures 

We now limit ourselves to the case v oo = and the two types of domains in Figure 
Eta)-(b), where <9Qi (the solid curve in (a) and the solid closed curve in (b)) is the boundary 
of the obstacle, 80,2 (dashed line segments in both (a) and (b)) is the far field boundary, 
and Q is the domain bounded by d£li and d£l2- This assumption implies 9 = on 80,2- 
In case (b), the circulation about the boundary is zero. 
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Figure 6. Domains 



Proposition 8.1. Let w e = (u £ ,v £ ) be a solution to (|5.1jl - (|5.2jl with Uoq > and Voo = 

on either domain £1 in Figure 6(a)-(b) ; satisfying that q £ < q* < q cav and 8 s is bounded. 
Then the integral 

ef^ (W)|W £ | 2 + CJ2 (p £ )^^|Vp e | 2 ) dxdy 
is bounded uniformly in all e > 0. 

Proof. We choose the special generator H* in (|7.7|) which yields V* of the form: 

v* = ^ + p( P ), v; = p'{ P ) = ^£ f q % v; = e. (8.i) 

^ H J q PQ. 

and (|7,1|) becomes 

r q ( pe ) da 

Ql(w £ ) x + Q* 2 (w e ) y = -e9 £ V(a 1 (p £ )V9 £ )+e J -± V(a 2 (p £ )Vp £ ), (8.2) 

where 

Q*i(w) = pqH* cos 9 — qHl smO, Q*2i w ) = pqH* sin 9 + qHg cos 9. 



Then 



Ql(w £ ) x + Q* 2 (w £ ) y = ediv^-a 1 (p £ )9 £ V9 £ + a 2 (p £ ){j_ 9 ^) Vp £ ) 



+ e* 1 (p £ )\V9 £ f-ea 2 (p £ ^ 2d ^ 



p £ q £ dp 6 

q£ dq 



ediv (-a 1 {p £ )9 £ V9 £ + a 2 (p £ )Vp £ j 



PQ 

! ^in)W9 £ ? + ea 2 {p £ )^Q^p £ ?. 

{p £ q e y 



1.3) 
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Integrating Q8.3JI over Q and using the divergence theorem, we obtain with q = > 
that 



(Q 



an 



l{w e ),Q* 2 (w £ ))-nds=e [ ( - a x ( P e )O e Ve e ■ n + a 2 (p £ )( T —)Vp £ -n)ds 

Jdn y Juoa PQ J 



+ e J L (p e )|Vr| 2 + a 2 (p £ )^-^\Vp £ \ 2 ^ dxdy, 
■=h + h, 

(8.4) 

where Q\,Q 2 depend only on (p £ ,9 £ ) and are independent of their derivatives. Thus, 
from the L°° bound, the left hand side of ()8.4[) is uniformly bounded for all e > 0. More 
specifically, from (J7.6|) and the formula of H* , we have the following: 

On <9f2i, (Q*, Q 2 ) • n = <f^ £ which is uniformly bounded in e > 0; 

On the horizontal part of <9ft 2 , £ = and \(Q*,Q 2 ) ■ n\ = \Q*\ = \q £ 6 £ \ = 0; 

On the vertical parts of dQ 2 , £ = and \(Ql,Q 2 ) ' n l = 1*3*1 = \p{uoo)uooH*(uoo)\ is 
uniformly bounded in e > 0. 

Using the boundary conditions (|5.2j) . one has 

I, = £ [ a 2 { f ^) Vp £ ■ rids + e [ a 2 { f °° ^) Vp £ -nds 
JdQi Juov PQ Jan 2 Juov PQ 

\p £ (u £ ,v £ )-n\{ f 1 ^-)ds, 
h Ju x PQ 



which implies that 1% is uniformly bounded due to the L°° bound, since the integrand is 
like glng near q = and is well behaved. Therefore, Q8.4JI implies that I 2 is uniformly 
bounded. □ 



In terms of a general generator H, equation (|7.1f) becomes 

Qi(w £ ) x + Q 2 (w £ ) y = eV(<7i(p e )V0 e )(p e i^ fl - H e )+eV(a 2 (p £ )Vp £ )(H^ + -H 96 ), 



or 



Qi(w% + Q 2 {w £ ) y = ediv(<ri(^)V^(^fl^ - H e ) + a 2 (p £ )Vp £ {H^ + -H e 

- ea^VO 6 • V(p £ H, ld - Hg) - ea 2 (p £ )Vp £ • V(i^ + 

^ (8-5) 
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The entropy pair (Qi, Q2) satisfies ()7.2(l which can be written in terms of H, 

dQi c 2 c 2 — 1 

— sin 9 {pH^o - H e ) - q cos 9 ^ — [ H » + -Hoe) , 



op pq Q P 

— — = qcos9 (pH^e - H e ) - pqsm.9 (H^ + -H ee ), 
oa p 

QQ2 c 2 c 2 — g 2 1 

— — = cos 9 {pH^e - H e ) - q sin 9 5 — (iJ„ + -H 9 e) ■ 

op pq q 1 v P 

q sin 9 (pH^e - H g ) + pq cos 9 (H^ + -Fee) • 



(8.6) 



Proposition 8.2. Assume that q £ (x,y) > ay on U <Z £1 for some constant ay > 0, in 
addition q £ (x,y) < q* < q cav ensured by the invariant regions. Then 

d x Qi±(w e ) + d y Q2±(w £ ) are confined in a compact set in 

for any entropy flux pair (Qi,Q2) generated by H £ C 3 . That is, hypothesis (A. 2) of the 
compensated compactness framework is satisfied for such entropy pair (Qi,Qz) generated 
byHGC 3 . 

Proof. For such an entropy pair (Qi, Q%) generated by H E C 3 , equation ()8.5|) is satisfied. 
By Proposition 18.11 we have 

Jf := ediv(a 1 (p £ )(p £ ^ - H e )V9 £ + a 2 (p £ )(H^ + ±-H e0 f)Vp £ 

in 7J _1 (i7) as e — > 0, and 

J| := -eai(p £ )Vr • V (///J^ - JT fl ) - ea 2 (p £ )V / 9 e • V(iI M + -^JT ee ) 

is in L l (U) uniformly in e. On the other hand, (C}i(w E ), Q2{w 6 )) are uniformly bounded 
which yield that Jf + Jf is bounded in W~ 1,co (U). Then Murat's lemma [31] implies that 
Jf + J| are confined in a compact subset of H~ l (U). □ 

As a corollary of Proposition 8.2, we conclude 

Proposition 8.3. Let the viscous viscosity fields w £ = (u £ ,v £ ) have the speed q £ (x,y) 
uniformly bounded in e > away from zero when (x, y) is away from the obstacle boundary 
dQi, i.e. there exists a positive e-independent a = a(5) — > as 5 — > such that q £ (x,y) > 
a(5) for any (x,y) G = {(x, y) E ^ : dist((x,y), dfli) > 5 > 0}. T/ien 

<9xQi±('W e ) + d y Q2±(w £ ) are confined in a compact set in flj~ (fi), 

/or any entropy pair (Qi, Q2) generated by H G C 3 . 

9. Convergence of the Vanishing Viscosity Solutions 

As noted earlier, the reduction of the support of the Young measure v is accomplished 
via application of the commutation identity (|6,2j) . Here we follow a technique of Morawetz 
|28| and use the entropy generators obtained via classical separation variables first given 
by Loewner [^Sj. Specifically, we look for H as the following two forms: 

H n (p,9) = F n (p)e ±md , 
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and 

H n (jji,9) = K n {p)e ±nd . 
Hence, from the generalized Tricomi equation (|7.5jl for H, we see 

9 M 2 -1 •■ 9 M 2 -1 

F n + n 2 F n = 0, K n -n 2 K n = 0, 

P P 

where = -j-. Substitution into the representation for Q\,Q2 as given in Proposition 

17.11 generates an infinite sequence of entropy pairs (Qi±,Q2±) associated with F n ; and a 
similar construction can be done for K n . 

First we apply the commutation identity to (Q^2,Q^±) associated with F n . Set 

i = igfl - QthQfl - Q ( 2 Y ) - (Q2 + - Q^')(Qi- - Q { il'). 

Then 



I = ( (F n pq cos 9 + inF n q sin 9)e m0 - {F'J q' cos 9' + inF' n q' sin 9')e 

x ((F n pq sin 9 + inF n q cos 9)e~ in(> - {Fjq'sm9' + inF'J cos 9' )e~ ine 
- ({Fnpq sm 9 - inF n q cos 9)e ind - {F'Jq 1 sin 9' - inF'J cos 9')e in6 '^j 
x ((F n W cos - inF n q sin #)e" irie - {F'Jq' cos 0' - mi^g' sin 9')e- ind ' 



With a tedious calculation, we obtain 

< v ® v', I > 

v®v\ 2inpq 2 (F n F n + F n F r , 



+ 2e in(e - Ql) (-F n pq cos 9F'J q sin 0' + n 2 F n q sin fli^g' cos 9' 
+F n pq sin 9F'J q cos 0' — n 2 F n q cos ^i^g' sin 0' 
+ 2e in(e - ' ) m (-F n/ og cos 9F' n q cos 0' - F n q sin fli^g' sin 0' 

—F n pq sin 9F' n q l sin 0' — F n q cos 9F'J q sin 6 1 ' ) ) , 



and then 



v <g> u', npq 2 F n F n + \F n F' n qq' sin(n(0 - 0')) sin(0 - 0')(Pp' + ™ 2 ) 
— nF n F' n pqq cos{n{9 — 9')) cos{9 — 9') 

v (g> v\ ^{qF n - q'F n )(pqF n - p'q'F' n ) 

1 o fn+l 
H — sm 



2 V 2 



0') J (^nF n F' n pqq' + F n F' n qq\pp' + n 2 )) 
I sin 2 (^(0 - 0')) (2nF n Fj q q' - F n F' n qq'{pp' + n 2 ; 
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Thus, from Q, 

< v ® u', -{qF n - q'F^)(pqF n - p' q' F' n 



+ \ sm 2 ~ 0')J {2nF n F' n pqq' + F n F' n qq'{pp> + n 2 ^ 

+ \ sm 2 (^(# " 0')) (2nF n F' n pqq' - F n F' n qq'(pp' + n 2 )) >= 0. 

Now interchange n and — n to get the n 2 -terms equal to zero. Hence 
< v ® i/, (qF n - q'F' n ){pqF n - p' q' F' n ) 

+ 2 ( sin 2 ( 2±±(6 - 9')) + sin 2 ( ^1(0 - 9')) ) F n F' n pqq' >= 0. ^ 



Then we follow verbatim from the argument of Morawetz |28j to lead the following 
convergence theorem. 

Theorem 9.1. Let i>oo = 0, \uoo\ < Qcav, 1 < 7 < 3. Assume that there exists a positive 
e -independent a(5) — > as (5 — > suc/i i/iai q £ (x,y) > a(S) for any (x,y) e!lj = {(^,y) £ 
SI : dist((x, y), dtli) > 5 > 0}. T/ien i/ie following hold: 

(i) T/ie support of the Young measure v Xi y strictly excludes the stagnation point q = 
and reduces to a point for a.e. (x,y) £ hence the Young measure is a Dirac 
mass; 

(ii) The sequence {u £ ,v £ ) has a subsequence converging strongly in L 2 oc (0) 2 to an en- 
tropy solution in the sense of Definition 7.1; 

(iii) The boundary condition (u, v) ■ n > on dtli is satisfied in the sense of normal 
trace in [S]. 

Proof. First, we choose F n ~ q~ n , with n sufficiently large and q < q cr , to show that the 
support of v cannot lie in both q < q cr and q > q cr as in Morawetz px] . If the support 
of v lies inside q < q cr , equation (|9.1j) or the argument of Chen-Dafermos-Slemrod-Wang 
applies to show that v is a Dirac mass. If the support of v is in q > q cr , Morawetz's 
reconstruction of DiPerna's argument in 27^ by using K n associated entropies (Qi±, Q2±) 
works to show again that v is a Dirac mass. The boundary condition on d$l\ follows the 
weak form of the second equation in (|2.1|) with R2 = eV(a2(p)Vp). The entropy inequality 
follows from (fg3|) . □ 

Remark 9.1. From the boundary condition Q5.2|) for the viscous problem, we find that, if 
one can show that eo"2(p e )Vp e • n — > in the weak sense on <9f2i as e — > 0, then we can 
conclude the desired boundary condition (u, v) • n = on dQ\ satisfied in the weak sense. 
This is the case when the limit solution has certain regularity along the boundary, i.e. any 
shock strength is zero at its intersection point with the boundary, as conjectured for the 
shock formed in the supersonic bubble near the obstacle (cf. Morawetz |26l [29 ) . 

Remark 9.2. For the first boundary value problem in Figure 6(a), if dO,± is sufficiently 
smooth, one expects that the transonic flow is strictly away from stagnation, which would 
expect that the viscous solutions stay uniformly away from cavitation in the whole domain 
Q. However, for the second boundary value problem in Figure 6(b), one expects that the 
transonic flow generally has a stagnation point on the boundary near the head of the 
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obstacle since the circulation about the boundary dVt 2 is zero. The assumption that there 
exists an e-independent positive a(6) — » as 5 — > such that q £ (x,y) > a(5), for any 
(x,y) 6 Q,5, is physically motivated and designed in order to fit both cases. 



10. Resolution of the Viscous Problem 

As we have seen in the previous sections, it is crucial that our viscous problem possesses 
smooth solutions w e = (u £ ,v £ ). We address the issue in this section. Consider the viscous 
problem: 

\v x -u y = eV(<7i(p)V6>), 



with 



(pu)x + (pv) y = eV(a 2 (p)Vp), ' 

01 = 1, 0- 2 (p) = 5 — for q > V2q cr , (10.2) 

and 01, 02 any smooth bounded, positive continuation for q < y/2q cr , and with the bound- 
ary conditions on the obstacle <9fh 

{V0-n = O on dQt, 
ea 2 Vp • n — \p (it, v) ■ n| = on t90i, (10.3) 
(«,«)- (lioc^oo) = on dQ 2 , with qoo < q cr , 

where p is a function of q given by Bernoulli's law (|2.2j) for 7 > 1 and H2.5|) for 7 = 1, and 
n is the unit normal pointing into the flow region on d£l. 
Introduce a new variable 

a{p) := [ P a 2 (0dZ. 



Since cr'{p) = (y 2 {p) > 0, we can always invert to write p = p(a), as well as q = q(a). The 
advantage of this change of dependent variable is that we now have the viscous terms to 
be eA6 and eAa. 

We use the polar velocity notation 

u = q(a) cos 6, v = q(a)sin0, 

and write the viscous system (jl0.1|) as 

j (q{a) sin 9) x - (q(a) cos 9) y = eA9, 
I (p(a)q(a) cos 9) x + (p(a)q(a) sin 9) y = eAa, 

which yields 



(10-4) 



J q'(a) sin 9 a x + q(a) cos 9 9 X — q'(cr) cos 9 a y + £7(0) sin 9 9 y = eA9, 
\ {p{ (J ) ( l{ (J )) 1 cos ® — p(^)(l(^) sin 99 x + (p(a)q(a))' sin y + p(a)q(a) cos 99 y = eAa, 

(10.5) 

with the boundary conditions on 90 1: 

fV0-n = O, 

\ [ 10 6l 

I eVcr • n — \p(a)q{a) (cos 6*, sin 6*) • n| = 0, 
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and the boundary conditions on dQ,2 ■ 

a = a OQ , 6 = 600, (10.7) 

where doc = cj(/9oo) and 0qc are constants. It is convenient to have the homogeneous 
boundary conditions on dO.2- 

a:=a-aoo, 9 := 9 - 9^; q(a) := q(a + (Too), p(ct) := + (Too). 

Hence, we have the following boundary value problem: 

'eA9 = q'(a) sm(9 + 9^) a x + q(a) cos(# + 9^) 9 X 

-q'(a) cos(# + 9^) a y + g(o-) sin(# + 9^) 9 y , 



eAa = (p(a)q(a))' cos(6> + 9^) a x - p(a)q(a) sm(9 + 9^) 9 X 



(lOi 



(10.9) 



'0„ in n, (10.10) 



+ (p(ct)<7(ct))' sin(# + #oo) o y + p(a)q(a) cos(# + 9^) 6 y , 
with the boundary conditions on dtli : 
ieVe-n = 0, 

[eVa ■ n - \p(a)q(a)(cos(9 + 9^), sm(9 + 0^)) ■ n| = 0, 

and (a, 9) = (0,0) on <9fi 2 - 

Consider the classical equation 

(9, a) = XT(9,a), < A < 1, 

which corresponds to the boundary value problem: 

sA9 = Xq'(a) sin((9 + 9^) a x + \q(a) cos(9 + 9^) 9 X 

— Xq'(a) cos(# + 9^) a y + \q(cf) sin(# + #00) 9 y in fi, 

eAa = \(p(o)q(o))' cos(9 + 6*00) a x - \p(a)q(a) sin(# + 9^) 9 X 

+ X(p(a)q(a))' sm(9 + 9^) a y + Xp(a)q(a) cos(# + 9^) 

eV9-n = on dQi, 

eVa ■ n = \\p(a)q(a)(cos(9 + 9 OQ ),sm(9 + #00)) ■ n| on dfl\, 

(a, 9) = (0,0) on <9fi 2 , 

which is the original system with e replaced by eA _1 . 

Now we recall some results from the linear elliptic theory. Consider the mixed Dirichlet- 
Neumann problem 

Aw = fi in Q, 

^ = f 2 on dto u (10.11) 
^w = f 3 on <9fi 2 , 

where is a bounded domain in M 2 with the boundary dfl = d^liUd^- We now introduce 
weighted Holder norms. For an arbitrary open set S and an arbitrary a > 0, we have the 
usual Holder norms || * ||a If <5 > and o + b > 0, we define 

!!j = {1 £ SI : dist(x, dQ.2) > 5} , 

and 



w 



1(6) 



sup 5' 



a+bi 



ir 



a,Qs ' 



21 



G.-Q. CHEN, M. SLEMROD, AND D. WANG 



and then denote C?L the set of all functions with \\w\\ ( a> < oo, and C*(fi) = C 2 (ftUd^i)n 
C(fi). 

Lemma 10.1 (Lieberman[2I], Theorem 2(b)). Suppose dili £ C 2+a , dQ,2 is nonempty, 
and Q satisfies the global H-wedge condition and the uniform interior and exterior cone 
condition on dtl. Then there is a unique solution to the linear elliptic problem (|1U.11|) and 

\H\ { Z<C(\\h\\^ + 11/21^ + 11/311,) 

for some (5 > and constant C > 0. 

Remark 10.1. Note that Lieberman's theorem is given for the boundary condition on d£l\: 
Mw := f3 t DiW + dw = /2, where d < 0. However, as noted by Lieberman on page 429 in 
|24| . the Fredholm alternative holds, and if dQ,2 is nonempty, the restriction d < may be 
relaxed to d < 0. The domains shown in Figure El satisfy the above conditions. 

We define T : (0, S) — > ((9, a) to be the solution map associated with solving the 
decoupled mixed Dirichlet-Neumann boundary value problem obtained from substitution 
of (9,S) into the right hand side of (|l().l()l) . For (0,2) E (C^)) 2 , the associated 
functions /i and /2 give finite values on the right-hand side of the Schauder estimates of 
Lieberman's theorem and fs = in our problem. Therefore, (9, a) stays in a bounded 
subset of (C^ a ) 2 which is a compact subset of (C( 2 ~+m) 2 ' and T is a compact map of the 

Banach space B = {C^^) 2 into itself. 

We recall one popular version of the Leray-Schauder fixed point theorem. 

Lemma 10.2 (Leray-Schauder Fixed Point Theorem). Let T be a compact mapping of a 
Banach space B into itself, and suppose that there exists a constant M such that \\W\\b < 
M for allW€B and all A S [0, 1] satisfying W = XTW . Then T has a fixed point. 

To apply the above theorem to our example, we first note that the only solution when 
A = is (#, a) = (0,0), hence it is certainly bounded in B. For A / 0, we study the 
solutions to our original viscous system with e = eX^ 1 > e. The invariant region argument 
we have given provides the uniform L°°(f2) bound independent of e. The gradient estimates 
from the entropy equality yield that y/iV(8, a) are in a bounded set of (L 2 (f2)) 2 , and hence 
a/eV (6, a) are in the same bounded set of (L 2 (Q)) 2 for all < A < 1. 

Now we examine the equations themselves: 

eA6 = A/i, eAa = A/ 2 , 

where /i, fi are in a bounded set of L 2 (Q). Hence, for < A < 1, eA(9, a) are in a bounded 
set of L 2 (0) and, by the C alder on- Zygmund inequality (Theorem 9.9 in Gilbarg-Trudinger 
EH), (6, a) £ W 2 ' 2 (n) . 

Next we differentiate the both sides of our equations with respect to x, y, and we get 
the equations of the form 

sA9 x = A/ia;, eA9 y = \f\ y . 

Since the right-hand sides are in a bounded set of L 2 (Q), so are the left-hand sides. We 
continue in this fashion to see that there is a constant M\ (independent of m and A) such 
that 

IKMiu + iivtMii^-w <Mi, 
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for all 1 < m < oo. 

Then Morrey's embedding theorem applied to our case yields V(0, a) bounded by con- 
stant M\ = M in (C S (Q)) 2 for any s > 1 and (9, a) bounded in B for any r > by the 
constant M\ = M, where the constant is independent of A. Thus, all possible solutions of 
(9, a) = XF(9,a), < A < 1, are bounded by constant Mi = M with M independent of 
A. Therefore, the hypotheses of the Leray-Schauder fixed point theorem are satisfied and 
the viscous problem has a solution, i.e. we have the following theorem. 

Theorem 10.1. The viscous problem has a solution (u £ ,v e ) for every e > such that 
<f <q* < Qcav f or 7 £ [1, 3) , where q* is a constant independent of e > 0. 

Acknowledgments. Gui-Qiang Chen's research was supported in part by the National 
Science Foundation under Grants DMS-0505473, DMS-0244473, and an Alexandre von 
Humboldt Foundation Fellowship. Marshall Slemrod's research was supported in part by 
the National Science Foundation grant DMS-0243722. Dehua Wang's research was sup- 
ported in part by the National Science Foundation grant DMS-0244487 and the Office of 
Naval Research grant N00014-01-1-0446. The authors would like to thank the other mem- 
bers of the NSF-FRG on multidimensional conservation laws: S. Canic, C. D. Dafermos, 
J. Hunter, T.-P. Liu, C.-W. Shu, and Y. Zheng for stimulating discussions. 



References 

[1] J. M. Ball, A version of the fundamental theorem for Young measures, Lecture Notes in Phys. 344, 

pp. 207-215, Springer: Berlin, 1989. 
[2] L. Bers, Results and conjectures in the mathematical theory of subsonic and transonic gas flows, 

Comm. Pure Appl. Math. 7 (1954), 79-104. 
[3] L. Bers, Existence and uniqueness of a subsonic flow past a given profile, Comm. Pure Appl. Math. 

7, (1954) 441-504. 

[4] L. Bers, Mathematical Aspects of Subsonic and Transonic Gas Dynamics, John Wiley & Sons, Inc.: 

New York; Chapman & Hall, Ltd.: London 1958. 
[5] G.-Q. Chen, Compactness Methods and Nonlinear Hyperbolic Conservation Laws, AMS/IP Stud. Adv. 

Math. 15, 33-75, AMS: Providence, RI, 2000. 
[6] G.-Q. Chen, Euler Equations and Related Hyperbolic Conservation Laws, In: Handbook of Differential 

Equations, Vol. 2, Eds. C. M. Dafermos and E. Feireisl, pp. 1-104, Elsevier Science B.V: Amsterdam, 

The Netherlands. 

[7] G.-Q. Chen, C. Dafermos, M. Slemrod, and D. Wang, On two-dimensional sonic-subsonic flow, Com- 

mun. Math. Phys. 2006 (to appear). 
[8] G.-Q. Chen and H. Frid, Divergence-measure fields and hyperbolic conservation laws, Arch. Rational 

Mech. Anal. 147 (1999), 89-118. 
[9] R. Courant and K. O. Friedrichs, Supersonic Flow and Shock Waves, Springer- Verlag: New York, 

1962. 

[10] C. M. Dafermos, Hyperbolic Conservation Laws in Continuum Physics, 2nd Ed., Springer- Verlag: 
Berlin, 2005. 

[11] R. J. DiPerna, Convergence of approximate solutions to conservation laws, Arch. Rational Mech. Anal. 
82 (1983), 27-70. 

[12] G.-C. Dong, Nonlinear Partial Differential Equations of Second Order, Transl. Math. Monographs, 

95, AMS: Providence, RI, 1991. 
[13] L. C. Evans, Weak Convergence Methods for Nonlinear Partial Differential Equations, CBMS-RCSM, 

74, AMS: Providence, RI, 1990. 



26 



G.-Q. CHEN, M. SLEMROD, AND D. WANG 



[14] R. Finn, On the flow of a perfect fluid through a polygonal nozzle, I, II, Proc. Nat. Acad. Sci. USA. 

40 (1954), 983-985, 985-987. 
[15] R. Finn, On a problem of type, with application to elliptic partial differential equations, J. Rational 

Mech. Anal. 3 (1954), 789-799. 
[16] R. Finn and D. Gilbarg, Asymptotic behavior and uniquenes of plane subsonic flows, Comm. Pure 

Appl. Math. 10 (1957), 23-63. 
[17] R. Finn and D. Gilbarg, Uniqueness and the force formulas for plane subsonic flows, Trans. Amer. 

Math. Soc. 88 (1958), 375-379. 
[18] D. Gilbarg, Comparison methods in the theory of subsonic flows, J. Rational Mech. Anal. 2 (1953), 

233-251. 

[19] D. Gilbarg and J. Serrin, Uniqueness of axially symmetric subsonic flow past a finite body, J. Rational 

Mech. Anal. 4 (1955), 169-175. 
[20] D. Gilbarg and M. Shiffman, On bodies achieving extreme values of the critical Mach number, I, J. 

Rational Mech. Anal. 3 (1954), 209-230. 
[21] D. Gilbarg, and N. S. Trudinger, Elliptic Partial Differential Equations of Second Order, Springer- 

Verlag: Berlin, 2001. 

[22] L.D. Landau and E.M. Lifshitz, Fluid Mechanics, 2nd Ed., Butterworth-Heinemann: Oxford, 1987. 
[23] P. D. Lax, Hyperbolic Systems of Conservation Laws and the Mathematical Theory of Shock Waves, 
SIAM: Philadelphia, 1973. 

[24] G. M. Lieberman, Mixed boundary value problems for elliptic and parabolic differential equations of 

second order, J. Math. Anal. Appl. 113 (1986), 422-440. 
[25] C. Loewner, Conservation laws in compressible fluid flow and associated mappings, J. Rational Mech. 

Anal. 2 (1953), 537-561. 

[26] C. S. Morawetz, The mathematical approach to the sonic barrier, Bull. Amer. Math. Soc. (N.S.) 6 
(1982), 127-145. 

[27] C. S. Morawetz, On a weak solution for a transonic flow problem, Comm. Pure Appl. Math. 38 (1985), 
797-818. 

[28] C. S. Morawetz, On steady transonic flow by compensated compactness, Methods Appl. Anal. 2 (1995), 
257-268. 

[29] C. S. Morawetz, Mixed equations and transonic flow, J. Hyper. Din". Eqns. 1 (2004), 1-26. 

[30] F. Murat, Compacite par compensation, Ann. Suola Norm. Pisa (4), 5 (1978), 489-507. 

[31] F. Murat, L'injection du cone positif de H^ 1 dans W^ 1 ' 9 est compacte pour tout q < 2, J. Math. 

Pures Appl. (9) 60 (1981), 309-322. 
[32] S. Osher, M. Hafez, and W. Whitlow, Entropy condition satisfying approximations for the full potential 

equation of transonic flow, Math. Comp. 44 (1985), 1-29. 
[33] Z. Rusak, Transonic flow around the leading edge of a thin airfoil with a parabolic nose, J. Fluid Mech. 

248 (1993), 1-26. 

[34] D. Serre, Systems of Conservation Laws, Vols. 1-2, Cambridge University Press: Cambridge, 1999, 
2000. 

[35] M. Shiffman, On the existence of subsonic flows of a compressible fluid, J. Rational Mech. Anal. 1 
(1952), 605-652. 

[36] L. Tartar, Compensated compactness and applications to partial differential equations. In, Nonlinear 
Analysis and Mechanics, Heriot-Watt Symposium IV, Res. Notes in Math. 39, pp. 136-212, Pitman: 
Boston-London, 1979. 

G.-Q. Chen, Department of Mathematics, Northwestern University, Evanston, IL 60208. 
E-mail address: gqchen@math.northwestern.edu 

M. Slemrod, Department of Mathematics, University of Wisconsin, Madison, Wl 53706. 
E-mail address: slemrod@math.wisc.edu 

D. Wang, Department of Mathematics, University of Pittsburgh, Pittsburgh, PA 15260. 
E-mail address: dwang@math.pitt.edu 



